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Abstract 

A new method for finding approximate solutions of the Ginsparg- 
Wilson equation is tested in 2-d. The Dirac operator is first constructed 
and then used in a dynamical simulation of the 2-fl.avor Schwinger 
model. We find a very small mass of the 7r-particle implying almost 
chirally symmetric fermions. The generalization of our method to 4-d 
is straightforward. 
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Introduction 

The last two years have seen a tremendous increase in our understanding of 
chiral fermions on the lattice (see [| for reviews of recent developments) . 
It has been understood that the anti-commutator characterizing the chiral 
symmetry in the continuum has to be replaced by the Ginsparg- Wilson 
equation § 

75 D + D 75 = D 75 D , (1) 

for the lattice Dirac operator D. Based on this relation, chiral symmetry can 
be defined on the lattice and at least for abelian fields chiral gauge theories 
can be formulated on the lattice ||. 

Bringing in the harvest from these developments based on the Ginsparg- 
Wilson equation now hinges on our ability to find solutions for this equation. 
So far two types of solutions are known: The perfect action approach to field 
theories on the lattice (see || for an extensive review) produces the fixed 
point Dirac operator H which is a solution of (p|). Although conceptually 
very beautiful - in addition to chirality on the lattice the fixed point oper- 
ator has other nice properties such as e.g. perfect scaling - constructing the 
perfect Dirac operator is technically very demanding and so far it has been 
implemented only in two dimensions [?]]-[§]. The other known solution to (|l|) 
is Neuberger's Dirac operator JTO] which goes back to the overlap approach 



to chiral fermions on the lattice [11|. Neuberger's construction provides 
a beautiful and surprisingly simple projection of any decent lattice Dirac 
operator onto a solution of the Ginsparg- Wilson equation. However, the nu- 
merical evaluation of the inverse square root necessary for Neuberger's pro- 
jection seems to be a challenging and computationally expensive problem. 
We should remark, that besides perfect actions and the overlap formulation 
also the domain wall approach which uses an additinal dimension gives rise 
to fermions with exact chiral symmetry |12| . 

Here we present a new approach to finding approximate solutions of the 
Ginsparg- Wilson equation, in particular we test the new method in 2-d. In 
4 dimensions the approach proceeds in exactly the same way and a detailed 
description of the new method in 4-d will be presented elsewhere fl3|| . 

The central idea is to expand the most general Dirac operator D on 
the lattice in a suitably chosen basis. Inserting this expression into the 
Ginsparg- Wilson equation renders a system of quadratic equations for the 
expansion coefficients. The basis for D can be chosen in such a way, that 
there is a natural cut-off for the expansion - the length of the paths allowed 
in the basis elements - and the remaining finite system of quadratic equa- 
tions can be solved numerically. In this letter we discuss the construction 
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of the above mentioned basis for two dimensions. We derive the system of 
quadratic equations and test the resulting Dirac operator in a dynamical 
simulation of the two-flavor Schwinger model. 

The Ginsparg- Wilson relation as a system of quadratic equations 

We begin our construction with denoting the most general 2-d Dirac operator 
on the lattice in the following form: 

3 

D*,v = E r « E c p Ii U i ( 2 ) 

<X=0 P&Px,y l ^P 

Here T a are the elements of the Clifford algebra, which in two dimensions 
has only 4 elements: Tq =1, Ti = 71, T2 = 72 and l"^ = 75. The 7^ are given 
by the Pauli matrices 71 = u\ , 72 = 02 an d 75 = 03 . To each generator T a 
and to each pair of points x,y on the lattice we assign a set of paths p 
connecting the two points. Each path is weighted with some complex weight 
Cp and the construction is made gauge invariant by including the ordered 
product of the gauge transporters U\ (g U(N) or SU(N)) for all links I of p. 
The action is then given by S = J2 x ,y ^x^x^ipy, where x and y run over all 
of the lattice. 

The next step is to impose on S the symmetries which we want to main- 
tain. Translation invariance requires the sets V"y to be invariant under 
simultaneous shifts of x and y and rotation invariance requires the terms 
to be invariant when rotating x and y with respect to each other. Parity 
((xi,X2) — > (xi, —X2)) implies that for each path p (with coefficient cp*) we 
must include the parity-reflected copy with coefficient Sp arity ■ c^ where the 
si g ns s parity are defined by J 2 ^al2 = sp* arity ■ T a . 

More interesting are the symmetries C and 75-hermiticityQ. It is easy to 
see that both of them imply a relation between the coefficient for a path p and 
the coefficient of the inverse path p^ 1 . Implementing both these symmetries 
firstly restricts all coefficients Cp to be either real or purely imaginary (for 
a = 0,1,2 they are real and Cp is purely imaginary). Secondly we find 
that the coefficient for a path p and the coefficient for its inverse p _1 are 
equal up to a sign s^ harge defined by CT a C = s" harge ■ T^, where T denotes 

1 75-hermiticity is defined by requiring 75D75 = for the iattice Dirac operator D. 
Together with C-invariance this property essentially requires that 71 and 72 come with 
derivative-like lattice terms, while 1 and 75 should come with scalars and pseudo-scalars 
respectively. Thus 75-hermiticity can be viewed as a lattice remnant of the properties 
necessary for proving the CPT theorem in Minkowski space. 
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transposition and the charge conjugation matrix in our representation is 
given by C = 72. 

When implementing all these symmetries we find that paths in our ansatz 
become grouped together in diagrams where, up to sign factors, all paths in 
a diagram come with the same coefficient. It is most convenient to represent 
the resulting Dirac operator using a pictorial notation. In Fig. [T] we display 
the leading diagrams in our expansion of the Dirac operator, i.e. terms with 
paths on single plaquettes and shorter than 2. The Dirac operator with this 
particular choice of terms will be referred to as D2 in the following. 




Figure 1: Paths and coefficients for the Dirac operator D2. 



We have furnished the Clifford algebra element 75 with an additional imag- 
inary i, such that the coefficients si,S2,sz,pi,vi,V2 are real. For the gen- 
erators 75,71 and 72 we assigned the relative signs for different paths in a 
diagram, according to the rules discussed above. These signs appear next 
to the paths. For the generator 1 all paths come with plus signs and we 
omitted these in Fig. 1. In order to be explicit about the interpretation of 
our pictorial notation we give as an example also the algebraic expression 
for the two 71-terms in D2: 

jU=±l 

+ v 2 si g n (^) [U^{x)U v (x + p) + U v (x)U^(x + i>)} 6 x+ p, + o iy . (3) 

It is important to remark, that the terms displayed in Fig. 1 are only the 
leading terms of an infinite series of diagrams. Diagrams with arbitrarily 
long paths have to be included, since it is known, that there are no ultra-local 
solutions of the Ginsparg- Wilson equation [14|. 
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The next step in the derivation is to insert the diagrammatic expansion 
of D into the Ginsparg- Wilson equation. To this purpose we define 

E = -D - 75D75 + 75^75^ , (4) 

and finding a solution of the Ginsparg- Wilson equation corresponds to hav- 
ing E = 0. The linear terms in (||) are easy to evaluate. The sandwich of 
D with 75 leaves the terms with I and ^75 in Fig. 1 unchanged, while the 
terms with 71 and 72 change their sign. These latter terms thus cancel when 
adding D and 75D75, while the former terms obtain a factor of 2. 

Finally we have to compute the quadratic part 75D75D. Here we obtain 
products of diagrams, each of them with their respective Clifford algebra 
elements. The multiplication of these terms proceeds in two steps (as can 
easily be seen by writing out an example using the algebraic expression 
(compare @) for our pictorial notation): In the first step the two elements 
of the Clifford algebra have to be multiplied. Since the Clifford algebra is 
closed the result is again a single element of the algebra. In the second step 
the paths of the diagrams under consideration have to be glued together. 
In particular one takes a path appearing in 75-D75 and continues it with a 
path appearing in D. When the latter path traces back some of the links 
of the former path, these links are removed from the diagram since the 
corresponding gauge transporters cancel each other. 

When now collecting all terms we obtain an expansion of the opera- 
tor E defined in (Q) in terms of our diagrams^]. The crucial step in our 
construction is to now view the diagrams in the expansion of E as basis 
elements (it is easy to see that they are linearly independent in an arbitrary 
background gauge field). Each of these basis elements is multiplied with a 
quadratic polynomial in the coefficients Si,Vi,pi. In order to have E = 0, 
and thus a solution of the Ginsparg- Wilson equation, we have to find the 
zeros of all these polynomials simultaneously, i.e. we have to solve a system 
of coupled quadratic equations. We have already remarked that one needs 
infinitely many diagrams (and coefficients Si,V{,pi) for an exact solution of 
the Ginsparg- Wilson equation and thus the discussed set of equations is an 
infinite set of quadratic equations for the infinitely many coefficients Si,Vi,pi. 
This infinite set of quadratic equations describes all possible solutions of the 

2 A short comment is in order here: Obviously £ is a hermitian operator, while the 
diagrams for 71 and 72 are anti-hermitian. These anti-hermitian terms, however, get 
transformed into their hermitian version when computing E along the lines sketched above. 
The hermitian versions of the 71- and 72-diagrams differ from the original terms by a 
different assignment for the relative signs of the paths. 
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Ginsparg- Wilson equation. Many interesting questions, such as how differ- 
ent solutions are connected, can be formulated in terms of this infinite set 
of equations. However, such questions will be pursued elsewhere |l3|]. Here 
we now concentrate on solving a subset of these equations and constructing 
approximate solutions of the Ginsparg- Wilson equation. 

Construction of approximate solutions 

From now on we will work with the truncated Dirac operator D2, repre- 
sented by the diagrams shown in Fig. 1. As already remarked above, such 
a truncation gives a meaningful approximation since, due to universality, 
we are only interested in Dirac operators where the coefficients of the paths 
decrease exponentially with the length of the paths. Of course this property 
has to be checked in the end. 

When deriving the expansion for E along the lines of the last section, 
now for a Dirac operator with only finitely many diagrams, we of course 
end up with only finitely many equations. For the case of D2 with its 6 
parameters s±, S2, ss,p\, v i, t>2 we find all together 45 quadratic equations. 
This system is overdetermined and we cannot satisfy all 45 equations with 
only 6 parameters. This property reflects the already mentioned fact, that 
there are no ultra- local solutions of the Ginsparg- Wilson equation If 
we were able to solve all 45 equations, E would vanish exactly, and we 
would have found a solution of the Ginsparg- Wilson equation with compact 
support. Thus, with our 6 parameters we can only construct an approximate 
solution of the Ginsparg- Wilson equation. For this purpose we will select a 
dominating subset of our 45 equations and solve only this subset. 

Before we come to this step, let's first discuss our system for trivial 
background field. In this case we can Fourier-transform the Dirac operator 
and the Fourier transform of any decent Dirac operator must obey D{p) = 
i p 1 + 0(p 2 ), i.e. the constant term must vanish and the linear term has 
to come with slope 1. Fourier-transforming our D2 and implementing the 
above condition gives the following two equations for the coefficients: 

= si + 4s 2 + 8s 3 , (5) 
= 2 Vl + 8v 2 - 1 • (6) 

These two equations must be implemented necessarily, such that we need 
only 4 of the equations from the expansion of E in order to fix the 6 param- 
eters s\, S2, S3,p\,vi,V2- The basic criterion for choosing the most dominant 
equations is clear: Equations corresponding to diagrams with shorter paths 
are more important, since the coefficients for diagrams with longer paths 
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become exponentially small. Thus we first choose the equations for the di- 
agrams with factors sx,s 2 and S3 in Fig. 1. Furthermore, it can be shown 
to all orders, that there is no diagram in the expansion of E which cor- 
responds to the v\ term of Fig. 1, since there exists no hermitian version 
of this diagram (compare the second footnote). Thus the remaining terms 
are all of length 2 and we could choose as our fourth equation either the 
equation for the p\ term or the equation for the hermitian version of the 
i>2 term. Here an interesting observation can be made: As we approach the 
continuum limit, the gauge transporters XJ\ approach 1. Thus the weights 
from the gauge transporters on the two paths leading to each point in the 
pi diagram of Fig. 1 approach each other as well. Since they come with 
opposite sign they cancel in the continuum limit. Thus this term takes care 
of itself in the continuum limit. The corresponding equation is not particu- 
larly important and we decided not to implement it at this order. The same 
argument can be repeated for the hermitian version of the v% diagram. We 
instead use as our fourth equation one of the equations corresponding to a 
length-3 diagram in the scalar sector. In particular we choose the term with 
a maximum number of different paths. This term is unique and we end up 
with the following set of equations, 

= s\ - 2si + As 2 2 + 8sj + 8pl + 4vf + 16v% , (7) 

= 2sis 2 - 2s 2 + 4s 2 S3 + <±viv 2 , (8) 

= 2 Sl s 3 -2s 3 + s 2 2 , (9) 

= S2S3 - viv 2 , (10) 

which together with (|5|) and (|6|) are used to determine the coefficients. It 
is straightforward to find solutions using some numerical solver. Although 
we found several solutions, only one of them provides a decent Dirac opera- 
tor. However, since the unphysical solutions are very far away in the space 
of coefficients, identifying the interesting solution was not a problem. The 
coefficients for our solution are given in Table 1. It is obvious that the coef- 



Sl 


S2 


S3 


Pi 


Vl 


V2 


1.26121 


-0.18471 


-0.06529 


-0.15745 


0.36940 


0.03264 



Table 1: Values of the coefficients in the Dirac operator D2. 

ficients decrease quickly in size as the length of the paths increases. In order 
to properly show exponential decay one would have to include more and 
longer diagrams. This would also allow the Dirac operator to even better 
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approximate a solution of the Ginsparg- Wilson equation. However, already 
with our relatively modest expense in additional terms we find very good 
approximate chirality. This property will be explored numerically in the 
next section. 



Results from the numerical simulation with D2 

We now choose our gauge group to be U(l). We first show two plots of the 
complex eigenvalues for our operator D2 in gauge field configurations from 
quenched ensembles at /? = 6 and = 4. 



(3 = 6 




P = 4 




Figure 2: Plots of spectra of D2 in the complex plane for quenched gauge 
field configurations at (3 = 6 (left-hand side) and (3 = 4. The circle indicates 
the spectrum for exact solutions of the Ginsparg-Wilson equation. 



It is easy to show, that (1) together with 75-hermiticity implies that the 
eigenvalues of an exact solution of the Ginsparg-Wilson equation lie on a 
circle of radius 1 around 1 in the complex plane. We find that the physical 
sector of the spectrum of D2 is very close to this limiting circle, while the 
doubler branches deviate considerably. This is not particularly disturbing, 
since the doublers decouple from the physical spectrum when driving the 
system to the continuum limit. Most importantly, however, we find that the 
physical branch intersects the real axis very close to the origin, rendering 
essentially massless fermions. This holds for both the (3 = 6 and the (3 = 4 
ensemble. We furthermore find, that also the fluctuations of the eigenvalues 
at the origin are much smaller than the fluctuations for e.g. the standard 
Wilson operator with a fine tuned mass term. This leads to error bars for 
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the 7r-mass (compare below) which are one order of magnitude smaller than 



with the standard Wilson action [15]. 

In order to further test the chiral properties of D2 we performed a sim- 
ulation of the massless 2-flavor Schwinger model with dynamical fermions. 
In the continuum this model can be solved explicitly and it contains an iso- 
triplet of massless 7r-particles and a massive iso-singlet r] (for the Schwinger 
model with flavor see e.g. Jl6|). When putting the model on the lattice a 
violation of the chiral symmetry will lead to a non- vanishing mass for the n. 
The smallness of m n thus is a measure for the quality of any approximation 
to the Ginsparg- Wilson equation. In our simulation we compute the masses 
and the dispersion relations for both the ir and the n. 

We compare our results to the numerical results for the two known ex- 
act solutions of the Ginsparg- Wilson equation, i.e. the overlap Dirac oper- 
ator and the perfect action. The data for these cases were taken from |J. 
Furthermore we compare our results also to the results for an alternative 
approximate solution of the Ginsparg- Wilson equation presented in jlT]] . In 
particular we compare with the truncated perfect operator augmented with 
an additional Clover term (TP+C). We remark, that also finite size effects 
can contribute to m^. However, we estimated these effects to be smaller 
than our error bars. Furthermore, all data presented in Table 2 were com- 
puted for the same lattice size (16 2 ) such that for our direct comparison of 
lattice results finite size effects are irrelevant. All data shown in Table 2 
have the same statistics (10 4 measurements). 





D2 


TP+C 


perfect 


overlap 


continuum 




0.0340(13) 


0.1304(27) 


0.0088(5) 


0.0038(23) 


0.0 


m v 


0.334(27) 


0.347(10) 


0.351(10) 


0.338(13) 


0.325730 



Table 2: it and rj mass at (3 = 6 for various lattice Dirac operators as well 
as their values in the continuum. 

In Table 2 we show, for (3 = 6, the masses for these different approaches 
and also give the continuum values. For the 7r-mass we find a value of 
0.0340 which is about four times as large as the mass obtained with the 
perfect action and about ten times as large as the mass from overlap. These 
Dirac operators are on the other hand much more expensive to simulate 
and still give a non- vanishing due to cut-offs or numerical errors. When 
comparing to the other approximate solution, the TP+C operator, we find 
that we have reduced the 7r-mass by a factor of 4, while the cost of the 
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numerical simulation remained the same. We do even better for the 77-mass 
where our data provides the best approximation of the continuum result. 

Let's now compare the dispersion relations. In Fig. |3| we show the tt 
and r\ dispersion relations of our operator D2 (filled circles) and compare 
it with the dispersion relations of the other lattice approaches (symbols) 
as well as with the continuum result (full curve). For the tt dispersion 




Figure 3: Dispersion relations for tt (left-hand side) and rj. The filled 
circles are the numerical data for our operator D2, the triangles represent 
the overlap operator, diamonds the perfect action and we use squares for 
TP+C. The symbols are connected with dashed lines to guide the eye. The 
full line is the continuum result. 

relation we find that the our operator leads to a nice linear behavior up to 
p ~ 1. For larger values of p the dispersion curve begins to deviate from 
the continuum value, although not as strongly as the overlap curve. Only 
the perfect action manages to obtain values close to the continuum result 
throughout. For larger momenta also the TP+C operator is relatively close 
to the continuum curve, however, for p —* the violations of chirality are 
obvious. 

For the 77 the picture is similar: Only the perfect action has a perfect 
dispersion relation. The other actions lead to dispersion relations which for 
increasing momenta deviate more or less from the continuum result. Again 
we find that our D2 curve lies between the continuum curve and the data 
for the overlap operator. 
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Summary 

In this letter we have explored a new approach to Ginsparg- Wilson fermions 
and tested the new method in 2-d. To keep the paper self consistent the con- 
struction of the Dirac operator was also given in 2-d, but the generalization 
to 4-d is straightforward and will be presented elsewhere ||13|| . 

The central idea is to construct an expansion of the Dirac operator in 
a suitable basis and to rewrite the Ginsparg- Wilson equation to a system 
of coupled quadratic equations for the coefficients. The basis elements are 
realized by diagrams of paths on the lattice, with the paths constituting a 
diagram being related by symmetries (rotations, C, P and 75-hermiticity) . 
Although the expansion in principle has to contain infinitely many terms, 
the method comes with a natural cutoff: Only Dirac operators where the 
coefficients decrease exponentially with increasing length of the paths are 
of interest. Thus we can work with a finite number of terms, rendering a 
finite system of coupled quadratic equations. The dominating equations of 
the system can be solved and used to construct approximate solutions of the 
Ginsparg- Wilson equation. 

We implemented this procedure and constructed the operator D2 which 
in addition to the terms appearing in the standard Wilson action only re- 
quires L-shaped terms of length 2. These terms imply only a moderate 
increase in the cost of a numerical simulation but drastically improve the 
chiral properties of the action. This was established in a numerical simula- 
tion of the 2-flavor Schwinger model. 

We are currently implementing the new approach in 4-d, where the con- 
struction proceeds essentially along the lines presented here. We expect 
that also in 4-d a considerable improvement of the chiral properties can be 
achieved with only a moderate increase of the cost of a numerical simulation. 

Acknowledgement: The authors thank Thomas Lippert and Uwe-Jens Wiese 
for insightful discussions and Christian Lang for important advice and checks 
of our equations. 
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